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SUMMARY 


The results of Initial calculations to determine the total amount of 
NO x deposited in the stratosphere by the Space Shuttle Solid Rocket Motors 
(SRM) are presented. Detailed calculations of the flow properties and chem- 
ical composition in the exhaust nozzle and plume were performed, accounting 
for such effects as gas/particle nonequilibrium and nonequilibrium chemistry 
in the nozzle, plume shocks (including the Mach disc) and nonequilibrium 
chemistry in the mixing/afterburning region and downstream of the Mach 
disc. 


The nozzle calculations show that, to within a factor of two, about 
4.5 lbm/sec of NO x leaves the two SRM's. The total amount of NO x depos- 
ited in the stratosphere is related to the amount leaving the nozzle via an 
Overall Plume Enhancement Factor (OPEF), whose value depends upon the 
influence of afterburning and shocks in enhancing the exit plane NO x mole 
fraction. Initial calculations show that the OPEF «2, indicating the mass 
flow of NO x in the plume to be » 10 lbm/sec at 30 km altitude with a possi- 
ble error factor of ? 4. The error bounds account for possible errors in 
input data (e.g. reaction rate coefficients), and the neglect of several phe- 
nomena (e.g. recirculating flow in the shuttle base region and intersecting 
plumes from the two SRM's) which may prove to be important. 

, : i 3 C; ' <•’ > 

For a vehicleOvelacity of .3750)ftS/:sec ,■ .therefore ,;;the NO x deposition 
rate in the stratosphere is about 2.7 x 10“ 3 Ibm/ft (4 g/meter), to within a 
factor of 4. 
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An accurate determination of the total amount of NO x deposited in the 
stratosphere by the Space Shuttle Solid-Rocket Motors (SRM) requires that de- 
tailed exhaust plume [NOx] ■< distributions be calculated, accounting for the 


effects of afterburning and plume shocks. These effects will increase the 
amount of NO x in the plume over its value at the nozzle exit plane, the amount 
of the increase being defined as the Overall Plume Enhancement Factor (OPEF). 
The total amount of NO x deposited in the .stratosphere? is, then, the amount 
"leaving" the exhaust plume. Figure 1 shows a schematic of the variation of 
Local Plume Enhancement Factor (LPEF)'i -as a function of distance from the 



FIG. 1 SCHEMATIC SHOWING DEFINITIONS OF LOCAL PLUME ENHANCEMENT FACTOR (LPEF) 

AND OVERALL PLUME ENHANCEMENT FACTOR (OPEF) 


* [NOx] is defined as a local concentration of NO x . 
t NO x consists of NO, NO z and N z O. 

t The LPEF is defined, as the mass flow of NO x at any axial station (m nq^) 

divided by the mass flow of NO x leaving the nozzle, (m NT ~.) ^NO 

- • co ‘ inu x ex. 

. is calculated from m NOx . « \ „ pu Y NO x r dr ’ where P> u and Y NO x are 

gas density, velocity and mas s .fraction of NO x , respectively and r is radial 
distance- from the axis. - 
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nozzle exit. . Note tha 


the~LPEF inc rea s e slo an a symptotic value , the OPEF , 
which defines the end en tire "NO x plume". There” are ho plume reactions which 
will decrease local values of [NO&]; plume reactions involving NO x far from 
the nozzle exit, ih<the cooler^regions oTthe* flow , convert NO to-NO z . 

The procedure followed in the present study has been to obtain initial 
estimates of [NO x ] at the nozzle exit plane and in the plume, via a series of 
simplified calculations which neglect; the effects of (i) two-dimensional flow 
and gas/particle interactions in the nozzle and (ii) shocks in the plume. These 
results give both an order.-of magnitude- estimate oT[jNO x ] in the nozzle and plume 
and the dependence of the predictions on uncertainties in some' of the; input data, 
e.g. reaction rate: coefficients and turbulent .'mixing rates. Next, a more de-; ? : 
tailed analysis was, performed in which the ;effects noted above wer.e'a.ccounteid 
for and the OPEF calculations made. This report presents the principal re- : 
suits of both the simplified and detailed calculations and conclude s'i$ith ...a first 
determination of the total amount of NO x deposited in the stratosphere .by the 
Space Shuttle Solid Rocket Motors. • .. . ' ' 


IX. INITIAL NOZZLE AND PLUME CALCULATIONS 


A. Nozzle 

The nozzle exit plane [NO x ] was initially determined via one-dimen- 
sional kinetic calculations. These calculations were made using nozzle 
pressure, temperature and velocity distributions calculated via a standard 
nozzle thermochemical equilibrium; code (from. AFRPL) as input to an Aero- 
Ghem-developed no n e qu Librium streamline .'codjea (NEQSLINE) . 1 Additional input 
to the code were the equilibrium composition in the combustion chamber 
(Table I) and a suitable chemical reaction mechanism and rate coefficients 
(Table II). The upper and lower bounds of the rate coefficients given in • _■ 

Table II are our best estimates based on the available experimental data. 

Figure 2 shows that the calculated NO x mole fraction* freezes near the -throat 
at a level about an order of magnitude greater than that predicted by the 
equilibrium calculation at the exit. (The 1-D EQUIL. curve was computed by 
the standard thermochemical equilibrium program.) The "STD. RATES (FULL)" 
curve was computed via NEQSLINE using all the reactions listed in Table II ' 
except Reactions D1-D6. A detailed examination of the reaction rates for 
individual reactions showed that reactions NQ2, NQ4, N05, NQ8, NO10, P3, 


* In the nozzle NO is the dominant constituent of NO^.; N 2 0 and N0 2 mole 
fractions are orders of magnitude less than that-of NO. 
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O 

DISTANCE FROM CHAMBER, ft 


r 

I FIG. 2 NC> MOLE FRACTION IN SPACE SHUTTLE SRM 

/ X 

I One-dimensional kinetic calculations assuming 

i gas/particle equilibrium 


*,T-- 
* v^V.- 
v-« 


P5, P6, P10 and Pll could probably be eliminated with little effect on [NO x ] . 
The "STD. RATES (REDUCED)" curve shows that, indeed, eliminating the above 
reactions does not significantly influence the results. Additional nozzle calcu- 
lations were made with the remaining NO x Production/Depletion and NO x 
Precursor reactions set to their estimated maximum and minimum rates re- 
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TABLE I. SPACE SHUTTLE SRM NOZZLE PROPERTIES 


Propellant Composition 

Wt. % 

AH° (kcal/mole) 

nh 4 cio 4 

69.6 

-70.69 

A1 

16.0 

0.0 

*“'6.884 ^10.089 O 0 . 2X8 No. 264 

12.0 

-12.0 

^6. 15^6. 97 Oj. 17 No. 03 

2.0 

-28.3 

Fe 2 03 

0.4 

-197.3 


Equilibrium 

Calculations 



Chamber 

Exit Plane 

Pressure , atm 

41.6 

1.08 

Temperature,. °K 

3400 

2330 

A/A* 


6.7 

Mass Fraction Al 2 0 3 y‘particle,s 

0. 30 

0.30 


Gas Composition, 

Mole Fraction 

CO 

2 . 49(- 1 ) 

co 2 

1 . 74(- 2) 

HC1 

1 . 43(- 1 ) 

H 

3 . 80(- 2) 

h 2 

2 . 7 7 ( - 1 ) 

OH 

9 . 1 1 (-3) 

h 2 o 

1 . 52(- 1 ) 

n 2 

9 . 9 2 (- 2 ) 

N 

6 . 19(-6) 

NO 

7 . 05(-4) 

O 

7 . 45(-,4) 

o 2 

1.62(-4) 

Cl 

1 . 30 (-2) 

HNO 

1 . 15(-6) 

n 2 o 

7 . 93(-8) 

NH 

2 . 78(-6) 

nh 2 

1 .07(-5) 
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TABLE n. SPACE SHUTTLE SR M NOZZLE/PLUME 
REACTION MECHANISM 


k f = A T~ N exp'(B/RT) 



i 

Heat Release 

A a 

N 

B 

(cal/ mole) 

Upper 

Error 

Bound 

Lower 

Error 

Bound 

HR1 

O +. O + M 


O z + M 

1 . 0(-29) 

1.0 

0 

30 

30 

HR 2 

O + H + M 

- 

OH + M 

1 ; 0(-29) 

1.0 

0 

30 

30 

HR 3 

H + H + M 

- 

H 2 + M 

2. 8(- 30) 

1.0 

0 

30 

30 

HR 4 

H + OH + M 

- 

H 2 0 + M 

2. 0(-28) 

1.0 

0 

10 

10 

HR 5 

CO + O + M 

- 

C0 2 + M 

1 . 0(-29) 

1.0 

-2500 

3 

30 

HR 6 

OH 4- OH 


h 2 o + O 

l.O(-ll) 

0 

-1000 

5 

5 

HR 7 

OH + H 2 

— ► 

H 2 0 + H 

5. 0(-l 1) 

0 

-5200 

3 

3 

HR 8 

O + Hfc\ 

- 

OH + H 

5 . 0(- 11) 

0 

-8460 

3 

3 

HR 9 

H + O z 

- 

OH + O 

3 . 0 (- 1 0) 

0 

-16500 

3 

3 

HR 10 

CO + OH 

- 

CO z + H 

5.0(-l 3) 

0 

-600 

3 

3 

HR 1 1 

H + Cl + M 

- 

HCl + M 

5 . 5(- 31 ) 

1 . 0 

0 

30 

30 

HR 12 

HC1 + OH 

- 

H z O + Cl 

7 . 2(- 1 2) 

0 

-3250 

30 

30 

HR1 3 

H + HC1 

- 

Cl + h 2 

8 . 8(- 1 1) 

0 

-4620 

10 

10 

HR 14 

OH + Cl 

- 

HCl + O 

3 . 0(- 1 1 ) 

0 

-5000 

30 

30 



NO Production/Depletion 






NOl 

o + n 2 


NO + N 

1. 3(-10) 

0 

-76000 

3 

3 

N02 

N + O z 

■«— 

NO + O 

2 . 2(- 1 1 ) 

0 

-6250 

3 

3 

N03 

NO + H ■ 


N + OH 

1 . 5(- 1 0) 

0 

-47000 

5 

5 

N04 

N + C0 2 

- 

NO + CO 

3. 0(-l 3) 

0 

-34000 

30 

10 

N05 

NO + NO 

- 

n 2 o + O 

2. 1 (- 1 2) 

0 

,-64000 

5 

5 

N06 

NO '+ H 2 

- 

HNO + H 

l.O(-H) 

0 

-54600 

10 

30 

NO 7 

HNO + Cl 

- 

HCl + NO 

3.0(-ll) 

0 

-1600 

10 

100 

N08 

NH + NO 

- 

N z O + H 

1 ; 0(- 1 2) 

0 

-4200 

100 

100 

N09 

NO + H + M 

- 

HNO + M 

5 . 5 ( - 2 9 ) 

1.00 

600 

3 

10 

NO10 

NO + M 

- 

O + N + M 

8. 3(-29) 

1.0 

-149000 

30 

30 


cm-molecule - sec units 

multiply A by upper error bound to get maximum possible value of rate 
coefficient 

divide A by lower error bound to get minimum possible value of rate „ 
coefficient 

(continued) 
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TABLE II 
(continued) 



NO Precursor 

A a 

N 

B 

(cal/mole) 

b 

Upper 

Error 

Bound 

Lower 

Error 

Bound 

PI 

N z O + H 


N 2 + OH 

1 . 3(-10) 

0 

-15200 

3 

10 

P2 

n 2 o + h 2 

•*— 

h 2 o + n 2 

1 . 0(- 1 3) 

0 

0 

100 

100 

P3 

n 2 o + CO 


n 2 + co 2 

3. 3(- 1 3) 

0 

-17400 

10 

10 

P4 

N z O + M 

■* 

N 2 + O + M 

6. 7 (-11) 

0 

0 

30 

30 

P5 

NH + Cl 

-* 

HC1 + N 

2 . 0 (- 1 1 ) 

0 

-4200 

10 

100 

P6 

NH + H 

- 

H 2 + N 

l.O(-ll) 

0 

-4200 

30 

100 

P7 

NH + OH 

- 

H 2 0 + N 

2. 3(- 1 0) 

0 

-2000 

3 

3 

P8 

NH 2 + Cl 

-► 

NH + HC1 

1.0(-10) 

0 

0 

3 

100 

P9 

nh 2 > h 

Z 

NH + H 2 

5 . 0 ( - 1 1 ) 

0 

0 

3 

100 

P10 

NH + H + M 

- 

NH 2 + M 

3. 0(- 30) 

1.0 

0 

100 

10 

Pll 

N + N + M 

t. 

N 2 + m 

5 . 6(- 30) 

1.0 

0 

10 

10 


D1 

no/no 7 

no 2 + no 2 - no+no+cv 

4. 0(- 1 2) 

0 

-27000 

3 

3 

D2 

NO + O + M 

- N0 2 + M 

3 . 0 ( - 3 3 ) 

0 

2000 

3 

3 

D3 

no 2 + O 

•z NO + O z 

1 . 7(- 11) 

0 

-600 

3 

3 

D4 

NO z + H 

- NO + OH 

5 . 8(- 1 0) 

0 

-1480 

3 

3 

D5 

NO + CO z 

- no 2 + CO 

3. 3(- 12) 

0 

-81600 

30 

30 

D6 

NO + O 

Z NO z + hfc 

6 . 6(- 1 2) 

2.0 

0 

3 

3 


- — M mm* - ...f- -■ 

j spectively. These results, also shown in Fig. 2»serve to give a preliminary 

estimate of the probable error bounds on the calculated [NO x ] at the nozzle 
j exit plane. 

. ' 

| B .' Afterburning Plume 


In order to get a ’feeling' for the potential increase of [NO x ] in the 
plume due to afterburning, a series of calculations were made with the Aero- 
Chem parallel mixing/chemistry (LAPP) c.dde 2 us ing . an. expanded plume 
boundary* as the dividing streamline. Initial conditions for the afterburning//- ?•> ' 
calculations (given in Table III) were obtained by assuming a chemically frozen 




The expanded plume boundary is defined here as f: -the radius x>f. the inv.iscid 
. boundary (dividing streamline) at which ambient pressure is Ve.ached.’' '' •. * 
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TABLE HI. INPUT DATA FOR SPACE SHUTTLE SRM PLUME 
AFTERBURNING CALCULATIONS 3 ALONG 
EXPANDED PLUME BOUNDARY AT 30 km 


Radius , ft 
Velocity, ft/sec 
Temperature, °K 
Pressure, atm 
Composition, Mole Fraction 
CO 

co 2 

HC1 

H 

h 2 

OH 

H z O 

n 2 

N 

NO 

no 2 

O 

o 2 

Cl 

A1 2 0 3 


Initial Conditions 


47.1 
1 .04(4) C 
1130. 

0.0118 


Free Stream 
Conditions 


3. 52(3) 
229- 
0.0118 


2 . 33 (- 1 ) 

2 . 09(-2) 
1. 54(-l) 

5 . 47(-3) 
2. 79(-l) 
4. 57(-4) 

1 . 41 (- 1 ) 

8 . 41 (-2) 
1.44 (-8) 

1 .44(-5) 

5 . 35 ( - 1 2) 
6 . 20 (- 6 ) 

1 . 32(-6) 

2 . 1 6 ( - 3 ) 
7. 99(-2) 


7. 90(-l) 


2 . 10 (- 1 ) 


Using the LAPP code 2 

Actual expanded radius is 33.4 ft. The radius used for these calculations 
is an 'equivalent' radius that accounts for the mass flow from both SRM 
nozzles . 

C A(B) = A x 10 B . 
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one -dimensional expansion. from the nozzle exit plane to the ambient pressure 
at 30 km, and solving simultaneously the continuity, momentum and energy 
equations. The initial' * equivalent' radius used .in th^plumeAcalctilations is 
nTJ larger than the calculated radius 'to account, for the total" mass flow from 
the two SRM ' s . 

The reaction mechanism used in the afterburning calculations includes . 
the HR1-HR14, NOl, N02, N03, Pll and D1 -D6 reactions of Table II. Typical 
results are shown in Fig. 3, which plots the ratio of the local [NO x ] to the frozen 
[NOx] (i. e. if afterburning were neglected) through out the -.plume Thi s demon- 



FIG. 3 INFLUENCE OF AFTERBURNING ON PLUME N0 x CONCENTRATIONS 


Space shuttle SRM; 30 km altitude; = 47. i ft; initial temperature = 1130°K 



y 


. t Note that .the Local Plume Enhancement Factor (LPEF) is related to the 
v >ar ea under the curve at. each axial station. (See definition on p. *1 ) ' ■ • • 

. v - * * rP . * 
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strates the direct effect of afterburning in enhancing plume [NOx] concentrations, 
the maximum increase being, a factorof about 2.4. If we arbitrarily increase 
the initial temperature from 1 1 30°K (used to obtaim - the results in Fig. 3) to 
2000 o K (to account for gross uncertainties ih this-simplified analysis)* 1 the , ‘ 

maximum increase in NO x over its frozen value is a. factor of about 35 ! • 'This 
result indicates the importance of accurately establishing the initial tempera- 
ture in this type of "expanded plume boundary" calculation. Additionalpara- 
metric .calculations were made to determine The effects, on [NO x ] , of uncer-, ■ . 
tainti.es in rate coefficients for reactions governing the direct production of 
[NO x ] and those controlling the radical concentrations! These results will 
be presented in the final report. 



C. Mach Disc 




The potential influence of the Mach disc on [NO x ] in' the plume was 
determined by using the NEQSLINE code, starting just downstream of'fhe Mach . 
disc with post- shock properties (assuming frozen chemistry Through the shock) 
as initial conditions . The. Mach disc location was estimated from the empirical 
..correlation of . Lewis and Carlson, 3 modified for free- stream velocity effects 
by DlAttore and Harshbarger , 4 For the purpose of- this-.,calculation the pressure 
distribution between the Mach disc and second plume shock was assumed to be 
similar to the pres sure distribution-f rom , the, nozzle exit to the throat. The.. ^ 
temperature was. then determined from a constant y calculation. 


Figur e ; 4 show s the results of the Mach disc calculation using both the . ; 
standard set of rate coefficients (Table II), and a reaction mechanism in which 
'the NO and P reaction rates (Table II) are set to their upper and lower: bounds. 
About an order of magnitude increase in [NO x ] is seen to occur as the flow 
passes through. the Mach disc (using the standard rates), due primarily to the 

h,gh temperatures (= 3000°K) in this region., . For reference. the chamb,r " 
value of lNO x J is also shown. This result indicates that the Mach disc will 
make an important contribution to local plume [NO x ] levels. 


III. OVERALL PLUME ENHANCEMENT 
FAG'TOR • G AL C U LA-TION S • 

The calculations discussed in Section II were directed at determining 
local effects and trends importaritih determining approximate levels of [NO x ] 
in the nozzle, and plume ./ Calculatioris aimed at more accurately determining 
the total amount of NOy in the plume are discussed -below. . 


Temperature was varied since [NO x ] is strongly dependent o,n iocal plume 
J; 'temperatures'. X* ’ - 
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distributions along the wall and centerline are given in Fig. 5, together with 
the one -dimensional- equilibrium results. These results are typical of bell- 
shaped nozzles, i.’:e. 'the flow near the wall rapidly expands just downstream of 
the throat,... giving lower wail pressures in this region., , The. exit plane wall 
pressure however is higher than the centerline pressure, i. e,.. p<^ = 7.2 psia 
and P wa ii = 21 psia. The gas temperature distributions. :do not' follow the pres- 
sure distributions because of the influence of solid particles. In fact, the 
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coupling -between' the gas and. particle- flows re, suits in the; exit plane gas tem- 
perature reaching a maximum valtie :(o'f 2 1 7 0 °K) of f - a xi s at rf r eX » 0 . 85 . 

Figure 6 com pare s,.:calculated values of NOx mole fractions' along the 
c enter line and wall, '('2- D kinetic)' with tho;se computed by the 1 -D - kinetic and 
1-D equilibrium programs'^, • The differ encef'be tween the 1-D kinetic and 2-D 
kinetic values at the exit ‘..plane is within the anticipated error bounds of [NO x ] 
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75 - 34 



B. Plume 

The procedure used to determine the enhancement of NO x in the plume 
due to afterburning and shocks was to (I) locate the plum e shocks via the 
AeroGhem unde r expanded rocket plume code (the AIPP code) 6 and (2) compute 
the NO x production using the LAPP code. 2 The shock structure shown in Fig. 
7 was scaled from a previous AIPP calculation of a Minuteman (MM) Stage 1 
plume in a similar flow region; The; LAPP code was then applied to three 
regions of importance; •(T) an -inner mixing zone (Region 1) in which the high 
temperature and pressure exhaust gas that has passed through the Mach disc 
mixes with lower temperature esdiaust products that have passed through the 

reflected shock; ( 2) an outer m ixing zone (Regio n 2) where cool exhaust 

products which have pa ssed ^through the exhaust: ga s shock mix with air that 
ha.s passed through the air shock; and (3) the merged shear layer (Region 3) 
which is initiated where (1) and (2) meet. Region 1 is characterized by an 
axial decay in pressure from the post- shock value to ambient. Regions 2 and 
3 are treated as constant pressure mixing zones. 
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FIG. 8 MERGED SHEAR LAYER LENGTH AND VELOCITY SCALES FOR USE WITH LAPP CODE 


1. Calculation Procedure 


The calculations for Regions 1 and 2 utilized the standard version of 

the LAPP code. Region 3, however, is initially characterized by two distinctly 

different shear Tayers.. The turbulent mixing within this region cannot therefore 

be adequately described by the single 'eddy viscosity normally provided for in the 

code. Consequently, ..the .program was modified to calculate separate character- 

istic length and velocity 'scales for the two shear layers until the position of the 
velocity maximum reaches the* centerline. (See Fig. 8.) At this point the pro- 
gram proceeds with the viscosity computation in -the normal manner. In addition, 
the empirical formulation for eddy viscosity developed by Stowell and Smoot, 1 * * * * * 7 
obtained by correlating data for reacting and non- reacting jets into moving sec- 
ondary streams, has been added to the code. Parametric runs have been made 
using both the Donaldson/Gray 8 and Smoot 7 eddy viscosity formulations. The 

reaction mechanism employed was the HR1-14, NOl-3, Pll set given in Table 
II. Thus far, NO — * N0 2 reactions have not been considered here, but will be 
treated in subsequent calculations. 
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TEMPERATURE, K 



Input data for the three regions are shown in Table IV. The input mole 
fractions for the Region 1 calculations were obtained by extending the non-, 
equilibrium' nozi le -c a 1 cula t io n s . .using to the To cation of the lvla'ch 

disc. The input for' Region 2 was scaled from the MM stage. 1 AIPP calcula- 
tion. The Region 1 calculation extended from ..the .Mach disc to^ambient pressure 
The extent of Region 2 was identical- to that of Region 1 The resulting profiles 
of temperature,, yplo city , and c ompo s ition were then used as input to the Region 
3 merged shear layer calculation. - , The velocity and -temperature profiles used 
as input data for the LAPP calculations are shown in Fig. 9 '. 

2. NO x Production Downstream of Mach Disc 

The total mass flow of NO x , at any ax'iaL, station was determined by radial, 
integration of the local values from the axis to the edge of the plume. The 
LPEF .was determined -by comparing the local value of [NO x ] with that of an 
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TABLE IV. INITIAL CONDITIONS FOR REGIONS 
1 , 2, AND 3 CALCULATIONS 


Region 1 


f 

Center - 

Free - 

line 

stream, 

Jet Radius, ft ! 17 


Vel. , ft/sec S 1250 

8800 

Temp., K 3370 

1200 

Press., atm 0.0608 a 



Mole Fraction! 


CO 

; 2.26C-1) 

1 . 76(- 1 ) 

co 2 

2 . 50(-2) 

3. 55(-2) 

HC1 

1 . 42(- 1 ) 

1 . 27(- 1 ) 

H 

7. 32(-3) 

l;59(-3) 

h 2 

2 . 76(- 1 ) 

2 . 24(- 1 ) 

OH 

6. 53(-8) 

2. 54(-6) 

h 2 o 

1 .45(-l) 

1 . 74(- 1 ) 

n 2 

9 . 38(- 2) 

2 . 60(- 1 ) 

N 

2 . 58(-9) 

1.0(-10) 

NO 

1 .61(-4) 

1 . 00(- 5) 

INERT 

1 . 6l(-4) 

1 .00(-5) 

O 

1 . 73(-8 ) 

1 . 97(-8) 

o 2 

5. 93(-7) 

1 . 49(- 7) 

Cl 

5 . 45(- 3) 

1 . 76(- 3) 

A1 2 0 3 

7 . 97(-2) ... 

. 0. 


Region 

2 

Region 3 

i 

C enter f 

F ree - 

C enter - 

1 ' Free- ' 

1 ine stream 

line 

stream 

36 


52 


9500 

3750 

5755 

3750 

1000 

220 

2233 

220 

0.0118 


0.0118 

| 


1 ♦ 76 (- 1 ) 


2. ll(-l) 

1 

3. 55(-2) 


2 . 1 2(- 2) 


1.27 ( - 1 ) 


9 . 97(- 2) 


1 . 59(- 3) 


l.ll(-l) 

i «* 

2 . 24(- 1 ) 


2,:25(-l) 

, 

2. 54(-6) 


8. 72( - 3) 


1 . 74(- 1 ) 


1 * 2 3 (- 1 ) 

, 7 . 9 ( - 1 ) 

2. 60(-l) 

7 • 9 ( - 1 ) 

8 . 70(- 2) 

1 .0{-10) 


2 . 05(- 6) 


1 .00(-5) 


3 . 0 1 ( - 4 ) 

1 

1 .00(-5) 


1 . 48(-4) 


1 . 97(-8 ) 


3 . 02(- 3) 

j 2 . 1 (— 1 ) 

1 .49(-7) 

2.1(-1) 

6 . 96(-4) 

1 . 76(- 3) 


3. 59(-2) 

j 

0. 


7. 36(-2) 

1 

1 

i 


Refers’ ‘to Tnitial:press;ure, flow is as.sume.d to expand in region 1 
to .ambient pressure 
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inert species having the same molecular weight and initial mass fraction profile. 
As mentioned previously, the overall plume enhancement factor, OPEF, is the 
value of the LPEF at the 1 end 1 of the NO x plume; i. e. , the point where the 
chemical production of NO x is negligibly small. 


The calculation to determine the initial effect of the Mach disc (Region 1) 
on the LPEF was made using the D/G viscosity model. Because of the uncer- 
tainty of applicability of viscosity .mod els in the merged she ar layer, parametric 
runs using both D/G and Smoot models were made to test for their effect on 
afterburning, and hence, NO x -production in Re gion 3. 


The OPEF results, shown in Table -V, indicate that the. mass flow of 
NO x leaving the plurhe.is: 30 -6 5 7° gr eater ithan that leaving the nozzle. Figure 
10 shows the axiM of the LPEF.^- It is apparent that’most of the 
NO x - production in the plume occurs imme dial e ly behind the Mach disc; little 
N O x production occurs in the afterburning region. Indeed, the amount of 


^ I T5-3J 



FIG. 10 LOCAL PLUME ENHANCEMENT FACTOR FOR NO 

x 
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TABLE V : . OVERALL PLUME ENHANCEMENT OF 


V is'c q s-ifey MoLd^ 4 - •;* V L V -* • 

— ..-v"v V*\ ; 

DonaMson/Gray 8 ' ^ . - *‘C. . 

(K .= 0 . 024 


QPEF • L^v,- 

' ■ ' -’r.fs ' • 

V i! 

i;;659 


Dona I'd son/ Gray 

(K-=:-0 .. 0 99 6 ) ; V , , 


'Smoot 9 ' Lf L 


■ r' 


ly632- 


tLSV- ^S’. r, A »V 




K = constant in eddy viscosity expression 



that the- high temperature behirid , the^Mach disc i s not maintained ; The dis- 
sociative- effects of. the' extreme nonequilibrium in thi s area tend to reduce ■ . 

.the temperature quickly thus inhibiting the formation of;. large amounts of 
NO^. Also, the mixing rate- has,' ; small- .effect on the OPEF. While the Smoot ; •' 


model mixes faster than the D/G model, the small amount of afterburning 
that does occur provides too low a temperature to cause marked NO x production. 


TV. 


TOTAL NO v DEFO5IT , IQM:1N^ITE;STRAT0SFHERE 


' : The maximum mass flow of NO x in the plume is the mass flow leaving -the 

the nozzle multiplied by the OPEF. The total mass flow leaving the two SRM 
nozzles is about 17 , 500 lbm/sec. With a mole fraction of NO x at the nozzle 
exit of 1.8(-4) 2), the mass flow of NO x at the exit is 4. 5 lbm/sec 2). 

The OPEF variation due to eddy viscosity models is relatively small (1.5-1/65). 
An additional factor of 2 variation is pos sible however because of uncertainties 
in calculating the Mach. disc size, *and because such factors as: (1) recirculating 
flow in the missile base region; (2) multiple Mach discs; (3) intersecting plumes 
were neglected in the present ^calculations.. ' CV. > 
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75 -33 



FIG. 11 INFLUENCE OF EDDY VISCOSITY MODEL ON PLUME 
CENTERLINE TEMPERATURES 


Our best estimate thus far for the NO„ mass flow in the plume at 30 km 
altitude is 10 lbm/ sec with a possible error of "-f 4. For a vehicle velocity of 
3750 -fps the NO x deposition is then 2.7 (-3) lbm/ ft (4 g/meter) (* 4). Our 
future efforts will be aimed at reducing the error limits on this result by 
examining the effects -of the above uncertainties (see Section V). 


V. PLANS 

In Phase II of the present study we plan to investigate several phenomenon 
that have been neglected in the calculations made to date, but which could signif- 
icantly effect the predicted ^amounts of NO x deposited in the stratosphere. This 
work has been divided into the following three tasks: 
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A. Interaction of Exhaust Plume with Flow Around Shuttle Body and in the 

Base R,ec ir..cnIatlon ;Re g io n . . .■/'/' , . , 

TJi.e peyfiprrned. thus fa r have .assumed that (a) me 

"free ^^^equal'to^the vehicle -Mach number' . 

.and- (b):.t^e^§ifeu^e .'base r eeiheulati on region is negligibly -small:.*- In; this task, 
we plan- to ..u|^2z;e>-e'xi.sting s'hufc^tle 'flow .field calculations at stratospheric 
altitudes tb d(feterni'ine actual flow properties approaching the base region. 

These properties w-ill be input to an analysis of the base recirculation region 
(e. g. Ref. 9)-which will be us J ed to determine inter alia the position of the , 
recompression shock and the composition of the recirculating base region 
flow entering .the air/rocket exhaust mixing region. Of prime concern here, is 
>• the possibility that’, because of the presence of a significantly large base re- ' • 
.circulation region, .the combustibles in the rocket exhaust products (H 2 and CQ) 
will ignite in the base region and, "feed", small quantities of; frde radicals (H, .. 1 ' 

OH) into the mixing region, thus considerably speeding up the chain- branching 
steps in the H 2 /O 2 reaction mechanism. The overall effect is that plume-'nfter- 
burni ng would be initi ated near the nozzle exit (i.e. the base would act as a 
■ flamehoider*) rather than far downstream where it would normally -occur due 
to mixing: between the exhaust products and ambient air .at "stratospheric alti- ■ 

tudes. The result would be- that .relatively large quantities of NO v could bd 
produced in the- plume- much-nealreir thernozzld-iexit ;plane than was predicted in 

Section III.- The base recirculation- region- propertied 1 will be input to the AIPP 
code and a calculation ^performed; to determine the total plume NO x production, 
including afterburning, and shocks . . This type of calculation should give ,a 
reasonably accurate quantitat-iye as.se ssment -of'the effect of the shuttle body/ 
base rec i mutation region on [NOx] 'levels -in-the plume . 

B. ' Intersecting-.^ . 

As' the shuttle altitude inc.r ea se-s;, the two solid propellant plumes become' 
greatly underexpanded and, in the stratospher'e/ t-hey will intersect close to the 
nozzle exit. The present plume calculations have been made by replacing the 
two intersecting plumes by a -single "equivalent" plxirne. having nozzle exit plane 
properties identical to those of a single nozzle, but. with the exit radius increased 
by V 2 to give the same total mass flow, t While this assumption is quite reason- 


* This "flameholding" effect has been observed in Apollo (first stage) flights. 10 

Apollo first stage ;( F - 1 engine), however , is a liquid rocket .using a LOX/ 

-.- RP-1 system. .V •’ '-d ■" . " 

t At lower altitude s : , e. g. about 10 km, it is probably. more reasonable to treat 
each plume separately. ___ 
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